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^ Abstract 

(— I ' In the beginning of the 1950's, Wigner introduced a fundamental deformation from the 

canonical quantum mechanical harmonic oscillator, which is nowadays sometimes called a 
Wigner quantum oscillator or a parabose oscillator. Also, in quantum mechanics the so- 
called Wigner distribution is considered to be the closest quantum analogue of the classical 
^ • probability distribution over the phase space. In this article, we consider which definition for 

such distribution function could be used in the case of non-canonical quantum mechanics. We 
then explicitly compute two different expressions for this distribution function for the case of 
I the parabose oscillator. Both expressions turn out to be multiple sums involving (generalized) 

^ ' Laguerre polynomials. Plots then show that the Wigner distribution function for the ground 

I state of the parabose oscillator is similar in behaviour to the Wigner distribution function of 

the first excited state of the canonical quantum oscillator. 
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^ . 1 Introduction 



The harmonic oscillator is one of the most appealing models due to its well known solutions 
in classical mechanics and in quantum mechanics, and due to a wide variety of applications in 
various branches of physics In classical mechanics its solution is unique, and so it is in the 
canonical non-relativistic treatment in quantum mechanics [2]. In this quantum approach, the 
position and momentum operators {p and q, respectively) satisfy the canonical commutation 
relation [p, q] = —ih, and the model is described by its Hamiltonian H and the Heisenberg 
equations. This algebraic description can be translated to a description by means of wave 
functions. In the position representation, for example, the operator q corresponds to multiplying 
functions by the variable q and the Heisenberg operator equations correspond to Schrodinger's 
equation, the solutions of which are expressed through Hermite polynomials. These solutions or 
wave functions lead to important position probability functions. 

To study the applicability of a system modelled by a harmonic oscillator, one often turns to 
expressions for the joint quasi-probability function of momentum and position, i.e. to a quantum 
analogue of phase space. The Wigner distribution function (d.f.) [3j is the main theoretical tool 
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here, as it is considered as the closest quantum analogue of a classical distribution function 
over the phase space. The quantum harmonic oscillator is one of the first systems for which 
the Wigner d.f. was computed analytically. Nowadays, there are also advanced experimental 
techniques that enable one to measure phase space densities for certain quantum systems [U [5] , 
and thus compare theoretical results with experimental results [6l[7l|8j. 

It is in this context that certain deformations or generalizations of the quantum harmonic 
oscillator are of relevance. In particular, deformations that are still exactly solvable and for 
which the Wigner d.f. can be computed analytically are highly interesting. One such model is 
the Q'-deformed quantum harmonic oscillator, being exactly solvable [HI \TU[ [TT] and having found 
applications in various physical models [12]. Its Wigner d.f. has been computed in [13]. 

A deformation or generalization of the canonical quantum harmonic oscillator that is more 
fundamental in nature was proposed by Wigner [14j. It is known as the one-dimensional Wigner 
quantum oscillator or as the parabose oscillator. In Wigner's approach, the canonical commuta- 
tion relations are not assumed. Instead of that, the compatibility of Hamilton's equations and 
the Heisenberg equations is required. This approach has been the fundamental starting point 
of so-called Wigner Quantum Systems [15]. In such an approach, one typically obtains several 
solutions (several algebraic representations), of which only one corresponds to the canonical 
solution. It is thus a natural and fundamental extension of canonical quantum mechanics. It is 
also of growing importance as deviations from the canonical commutation relations (in partic- 
ular of commuting position operators in the context of non-commutative quantum mechanics) 
are supported by investigations from several field [TBI [13 O US] . 

The one-dimensional Wigner quantum oscillator itself was already solved by Wigner |14j . 
leading to an oscillator energy spectrum similar to the canonical one but with shifted ground 
state energy equal to some arbitrary positive number a. In a sense, this number a can be 
considered as the deformation parameter, and the solutions of the Wigner quantum oscillator 
reduce to the canonical one when a = 1/2. Wigner's formulation was also the inspiration for 
the introduction of parabosons (and parafermions) , as the common boson commutator relation 
of the canonical case is, in Wigner's approach, replaced by a triple relation that is also the 
defining relation of a paraboson ^20j. Later on, one realized that Lie superalgebras are the 
natural framework for this algebraic description [21], and in particular that the representations 
of the Wigner quantum oscillator are in one-to-one correspondence with the unitary irreducible 
lowest weight representations of the Lie superalgebra 05p(l|2). 

So the Wigner quantum oscillator, being the generalization of the canonical one, could play 
an important role to explain various phenomena more precisely, and its importance is undeniable. 
However, so far its application has been restricted in comparison with the canonical harmonic 
oscillator due to the fact that its behaviour in phase space has not been studied. In other 
words, there is no expression of the joint quasi-probability distribution function (Wigner d.f.) of 
position and momentum available for the Wigner quantum oscillator. The main problem is that 
the treatment a Wigner d.f. usually starts from the assumption that the commutator of p and 
g is a constant, which is no longer the case for the non-canonical Wigner quantum oscillator. It 
is precisely this problem which is being tackled and solved here. 

In this paper, we first recall some basic facts of the one-dimensional Wigner quantum oscil- 
lator or parabose oscillator. The emphasis is on the algebraic description, in terms of the Lie 
superalgebra osp(l|2) and representations characterized by a positive real number a (section 2). 
In section 3 we discuss the main definition of the Wigner d.f. in the current non-canonical set- 
ting. The main part of the paper is devoted to computing the Wigner d.f. Wn{p, q) analytically, 
for the parabose oscillator in the nth excited pure state \n). For n = (and n = 1), the calcu- 
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lation is somewhat simpler than in the general case, and this is treated first in section 4. The 
calculation involves various computational tools. In some steps, series are term wise integrated, 
so one should pay attention to convergence. Here, we shall often assume that the parameter a 
is of the form a = 1/2 + m, with m a nonnegative integer; in that case all series appearing are 
finite and there are no convergence problems. The final expression of WQ{p,q) is essentially a 
series in terms of Laguerre polynomials, depending only upon p'^ + q^, or alternatively a conflu- 
ent hypergeometric series in + q"^, see (|33p . Section 5 then deals with the general case: the 
computation of Wn{p,q)- Here, the analysis is quite involved, and some intermediate calcula- 
tions are performed in an Appendix B. We obtain two alternative expressions for Wn{p, q), (|39p 
and ()4ip . both as series of Laguerre polynomials in p^ + q^ . Section 6 deals with some plots 
of the newly obtained Wigner distribution functions, discusses some properties and summarizes 
some conclusions. 

To end the introduction, let us introduce some of the notation and classical results we are 
going to use. Many of the formulas encountered in this article will involve binomial coefficients, 
factorials etc. It is therefor convenient to use the notation for (generalized) hypergeometric 
series |22] . Let r and s be nonnegative integers, a hypergeometric series with r numerator and 
s denominator parameters is then defined as 

I ai, . . . ,ar \ _sr^ (ai, . . . ,ar)k 

where 

Y(a + k) 

(ai, . . . , ar)k = (ai)fc(a2)fc • • • (ar)fc, and {a)k = . = a{a + 1) ■ ■ ■ {a + k - I) (2) 

r(a) 

is the rising factorial or Pochhammer symbol. When one of the numerator parameters aj is a 
negative integer, the series is terminating and in fact a polynomial in z. Of course, if one of the 
denominator parameters is also a negative integer, it has to be "more negative" such that the 
numerator vanishes before the denominator does. 

All classical orthogonal polynomials can be written using this hypergeometric notation. In 
particular one has for the Laguerre polynomials |23j : 

K^iz) = j iFi ,,;z]= j iFi{-n;a + l;z). 3) 

Here, we also introduced an alternative notation for the confluent hypergeometric series iFi. If 
the parameter a is 0, we will simply write Ln{z) instead of Ln\z). Often, the Ln{z) are called 
Laguerre polynomials, while the L^"'* are called the generalized Laguerre polynomials. For the 
classical orthogonality of the Laguerre polynomials to hold a needs to be bigger than —1. In 
particular, if a is a negative integer, the deflnition ^ is no longer valid since one has a negative 
integer appearing as a denominator parameter. The Laguerre polynomials can be written in an 
alternative way such that they are deflned for all values of a: 

Li-\z) = 1 V tlM^±±±lhzl,K (4) 

Many well known summation formulas can be very neatly expressed using the hypergeometric 
notation. Most notable are the binomial summation theorem: 

iFoia- -z) = (5) 
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with \z\ < 1 if a is not a negative integer, and the Chu-Vandermonde summation theorem for a 
terminating Gauss hyper geometric series: 

(-n,a \ {c-a)n , . 

with n a nonnegative integer. There exists also a very extensive theory of hypergeometric 
transformations, a very important transformation being Kummer's transformation for confluent 
hypergeometric series: 

iFi(6;c;z) = e^Fi(c-5;c;-z). (7) 



2 The one-dimensional parabose oscillator 

In this paragraph we briefly describe the treatment of a harmonic oscillator as a Wigner quantum 
system. Although first introduced by Wigner almost 60 years ago, and described in some articles 
and book chapters, it may still be useful to recall the main results. 

Consider the Hamiltonian for a one-dimensional harmonic oscillator (in units with mass and 
frequency both equal to 1): 

H = ^- + ^-, (8) 

where p and q denote respectively the momentum and position operator of the system. Wigner 
|14j already noted that there are other solutions besides the canonical one if one only requires 
the compatibility between the Hamilton and the Heisenberg equations. Thus, one drops the 
canonical commutation relation [p, g] = —i (we have taken ^ = 1), and instead one imposes the 
equivalence of Hamilton's equations 

A _ dH x_dH 
^ dq ^ ^ dp 

and the Heisenberg equations 

p = i[H,p], q = i[H,q]. 
Expressing this equivalence leads to the following compatibility conditions: 

[H,p]=iq, [H,q]=-ip. (9) 

So, one has to find operators p and q, acting in some Hilbert space, such that the compatibility 
conditions ^ hold, with H given by ([5]). Also, in this Hilbert space p and q have to be self- 
adjoint. 

The solutions to Q can be found by introducing two new (yet unknown) operators 6+ and 

b^ = ^{q^ip), (10) 

or equivalently 

q=l=ib^ + b~), p=±={l+-b-). 
It is then easily checked that 

H=\{b-M], (11) 
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and that the compatibihty conditions ^ are equivalent with 

[{r,6+},6±] = ±26±. (12) 

The relations (|12p are in fact the defining relations of one pair of parabose operators [20] . 
Moreover, from the self-adjointness of the position and momentum operators it follows that 

{b^)^ = y^. (13) 

It is known that the Lie superalgebra generated by two odd elements subject to the restric- 
tion ([12]) is the Lie superalgebra osp(l|2) [21j. So, the solutions to our problem are given by the 
star representations of the Lie superalgebra osp(l|2). These are known, and are characterized 
by a positive real number a and a vacuum vector |0), such that 

6-|0) = 0, {6-,fe+}|0) = 2a|0). 

The representation space can then be shown to be the Hilbert space £^(Z+) with orthonormal 
basis vectors |n) (n G Z+) and with the following actions: 



(14) 



S+|2n) = v^2(n + a) \2n + 1), b~\2n) = V2n\2n - I) , 

6+|2n + 1) = v^2(n + l)|2n + 2), 6"|2n + 1) = ^/2{n + a) |2n), 

from which it immediately follows that 

{6",6+}|n) = 2(n + a) |n). (15) 

Using the actions ()14p . one can now compute the action of the commutator [p,q\: 

[p,q]\2n) = -2ai\2n), [p, g]|2n + 1) = -2(1 - a)i|2n + 1). (16) 

So the canonical commutation relations are only satisfied for a = 1/2 and any representation 
with a 7^ 1/2 thus leads to a non-canonical solution. It is also important to note that for such 
non-canonical solutions, the commutator [p, q\ is not constant. This last fact has a serious impact 
on the computation of exponential operators involving linear combinations of p and q since the 
Baker-Campbell-Hausdorff theorem does not give rise to a simple formula. 

The energy spectrum is now very easy to determine as well using (I15p and (llip : 

H\n) = {n + a) |n), 

so one has an equidistant energy spectrum with ground level given by a. This once again confirms 
that only a = 1/2 yields the canonical solution. 

Also the wave functions for the parabose oscillator in the position and momentum repre- 
sentation are known. Working in the position representation, the operator q is still represented 
by "multiplication by g" , the operator p however is no longer given by its canonical realization 
—i-^, but instead has a different realization with an extra term (depending on a) coming in. 
Using this realization the time-independent Schrodinger equation can then be solved yielding the 
following expressions for the orthonormal wave functions in terms of Laguerre polynomials [241 
Chapter 23] and [25]: 



(17) 
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Again, for a = 1/2, one recovers the very well known expression for the wave functions for a 
harmonic oscillator in terms of Hermite polynomials, since there are formulas relating L^^^'^\x'^) 
to H2n{x) and H2n+i{x) respectively. 

3 Distribution function for a non-canonical quantum system 

A major advantage of the phase-space formulation of quantum mechanics is the fact than one 
is no longer dealing with operators, but instead with constant-number equations. More in 
particular, the expectation value of an arbitrary operator A{p, q) can be calculated using the 
distribution function (?) as follows: 

TY{p{p,q)A{p,q)} = jj A{p,q)^{q,p)dqdp, (18) 

with A{p, q) the scalar function obtained by replacing the operators p and q in the expression for 
A^p, q) by scalar variables p and q. Herein, p is the density operator, representing the state of the 
system. Also, both integrations run from — oo to -|-oo; this will be the case for all integrations 
in this article. 

It is also known that the distribution function <I> depends on the so-called correspondence 
rule used, this is, on the way one associates a function of non-commuting operators (thus an 
operator) to a given function of scalar variables. A very popular association rule is the Weyl 
correspondence, yielding as distribution function the so-called Wigner distribution function. 
Very often, the definition for the Wigner d.f. is written as follows: 

W{p.q) = ^ l{q+P\q-^)e~^^Pdrj. (19) 

This definition, however, is completely unsuitable for the case of non-canonical quantum me- 
chanics, since it relies heavily on the canonical commutation relations. Indeed, it was shown 
by Tatarskii [26] (see also section 2.1 of Lee's review paper [27J) that the Wigner d.f. defined 
by first principles (see equation ([20|) here) gives rise to expression (fT9|) only if p and q satisfy 
[j5, q] = —i. So if the canonical commutation relations are not satisfied, one should go back to 
use this first principles approach, which is described in the rest of this section. 

Another approach is the following [26]: is it clear that in order to define an association 
rule between functions of scalar variables f{p,q) and operators f{p,q) it suffices to define this 
correspondence for the function F{X,fj,) = exp{i{Xp + pq)), with A and p real variables. Indeed, 
an arbitrary analytical function f{p,q) can be obtained from F{X,iJ,) in the following way: 

Note that this expression is well-defined since the partial derivate operators commute. The 

operator corresponding to the function f{p,q) under the chosen association rule, i.e. for the 
particular choice of F, is then given by 

I d 1 d - 

f{p,q) = f{--^,--^)FiX,fi)\x=i,=o ■ 
I oX I op 

Once a choice for F is made, the corresponding distribution function Wp may be calculated 
as follows: 

WF{p,q) = ^ll {^mX,pM)e~-'<^'P+^^''UXdp, (20) 
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where ip denotes the state of the system. This definition is independent of the canonical com- 
mutation relations and still allows the computation of averages by use of (fTSj) . 
We now make the following choice for F 126\ [27] : 



(21) 



where we have introduced the operator X: 



X = Xp + fiq. 



(22) 



Since this is the only correspondence rule we are using in this article, the subscript F will 
be dropped from the notation. On the other hand, as we will be concentrating on calculating 
explicit expressions for the distribution function W for the pure states \n), we will introduce the 
subscript n to indicate this: 



Since in this article our goal will be to compute an explicit expression for (|23p , it is convenient 
to recall the expression for the Wigner d.f. Wn in the case of canonical commutation relations: 



where Ln are the Laguerre polynomials. So, whatever result we obtain for Wn, replacing the 
parameter a by 1/2 should yield (|24|) . 

4 Wigner distribution function for the ground state and the 
first excited state 

In this section, we compute an explicit expression for Wo{p,q). This will already show the 
techniques involved in the computation of the Wigner d.f. for arbitrary states |n), but the 
resulting calculations are not quite so lengthy. To reiterate, the goal is to compute 




(23) 




(24) 



vr 




(25) 



Thus, obviously, in a first step we are going to determine the matrix elements (0| exp(iX)|0) 
and in a second step the integration will be performed. 
The operator ex.p{iX) being defined as 




fc>0 



we first try to find (OlX'^IO). To this end, we write X = Xp + fiq in terms of b^: 

i> +?+ -7"- -1 + /U ± 
X = a b + a with a = ^— . 




(26) 



It is then immediately clear that 



(0|X 



-2k+l 



|0)=0 



(27) 
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since an odd number of applications of 6"*" and b" on |0) does not yield a term corresponding 
to |0) in the resulting expansion, as can be seen from (jl4p . On the other hand, using (jl4p . one 
finds easily that 

{o\x'\o) = 1, 

(0|X2|0) = 2a+a-a, 
(0|X^|0) = {2a+a-fa{a + l), 
{0\X^\0) = {2a+a-fa{a + l)(a + 2). 

This suggests, as will be proved later, that one has 

(0|X2^|0) = (2a+a-)'= (a)^ = (A^ + ^2)^ (28) 

Using (HZ]) and (I28|), one thus finds that 



fc! ' ' ^ (2A;)! ' ' ' ' ^ (2A;)! 

fc>0 A:>0 ^ ^ fc>0 ^ ^ 

(«)fc . + /a A^ + /i2 \ 

where we have used the fact that (2A;)! = (1/2)^ kl4^. This means that we now have: 



(29) 



47r2 yy ' ' vi/2' 

In order to compute this integral, it is useful to introduce a quadratic decaying exponential 
into the integrand by applying Kummer's transformation (j7]) on the confluent hypergeometric 
series. One then has: 

where we have swapped summation and integration. We assume that the parameter a is such 
that this is possible. In any case, when the parameter a is of the form 

a = 1/2 + m with m G Z+, (30) 

this can certainly be done since one is then only dealing with finite summations. In Appendix 
A, it is shown that 

4^// ^"^^''"''^^'^"'^^''^'''^(^' + /^')'^^^^ = ^e~^'~'''A;!4'=L,(p2 + (31) 

such that 

fc>0 ^ ' 
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It is important to remark that, apart from the exponential factor, the expression (|32p is a 
polynomial in + whenever the parameter a is of the form (j30p . In this case, it is possible 
to simplify the expression (j32p even further. Expand the Laguerre polynomials as single sums 
and switch the two (finite) summations. On the resulting inner sum Vandermonde's summation 
theorem ([6]) can be applied. After some further easy manipulations one then finds: 

W^{p,q) = le-P'"^'^^ + ^2)^ ^-^^^ a =1/2 + m. (33) 

Clearly, in the canonical case (i.e. m = 0) the resulting expression reduces correctly to WQ{p,q). 

To compute Wi{p,q), we first have to find an expression for the matrix elements (llX'^ll). 
Again, one finds the matrix elements with odd powers of X to be zero. For the even powers, 
one finds using and (p6]) : 



(1|X0|1) = 1, 

= 2a+a-(a + 1), 

= {2a+a-f{a + l)(a + 2), 

{l\X^\l) = {2a+a-f{a + l)(a + 2)(a + 3), 



indicating that 



(i|x^|i) = (o|x'=|o)U^, 



+1- 



Since these matrix elements are the only place where the parameter a appears, the expression 
for Wi{p^ q) is also found by replacing a with a + 1 in the expression for Wq{p, q): 

Wi{p,q) = Wo{p,q)\a-,a+l. (34) 



5 Wigner distribution function for arbitrary states \n) 

We now would like to find an explicit expression for the following integral: 

As is clear from the previous section, the key in calculating this integral is finding an expression 
for the matrix elements {n\X^\n). As before, matrix elements with odd powers of X are easily 
seen to be zero. It will also be the case that matrix elements for a state |2n+l) will be determined 
from those with state |2n) simply by replacing a with a + 1. We will now give two different 
expressions for the matrix elements {2n\X'^^\2n). We will show that they are equivalent and will 
compute the Wigner d.f. using each of them. 

The proof that the proposed expressions are in fact the correct matrix elements is relegated 
to an appendix, despite the importance of the result to the calculations. The proof, although in 
itself not difficult, requires rather lengthy calculations as one also needs the off-diagonal matrix 
elements. Basically, it involves checking that a certain expression satisfies the recurrence relation 
for the matrix elements (as well as the correct boundary conditions). 

We have the following expressions for the matrix elements {2n\X'^^\2n): 

(2n|12fc|2n) = (A^ + /.2)fc ^ M |^^^ (« + „ _ j)^^a + 2n)k-2j, (35) 

2j<k 



9 



or alternatively 



min(n,fc) ^ / \ 

(2n|12fc|2n) = (A2 + ^y ^ _ + j)k^,{k + 1 - j^j. (36) 

Note that both expressions reduce to (128p for n = 0. To see that both expressions are indeed 
equivalent, it is easiest to write them in hypergeometric notation. For formula (j36p we have 

{2n\X''\2n) = (A^ + f,')' {a)k 3^2 f ^+ ^' ""^ ; 1 



1, a 



^2fcio^\ m2 , ,,2Nfc/ , 17 -n, -A:/2, 1/2 - A:/2, 1 - n - a 



whilst formula (l35l) is 



(2n|X-|2„) = (A' + M')' (a + 2„)» j_ _ _ 1 1 „ - <./2 - ' 

So, we would like a transformation for hypergeometric series that relates a series of type 3-F2 
with one of type 4-F3. In Krattenthaler's list [28], transformation T3231 does just that and is 
indeed applicable to the 3-F2-series at hand. It will turn out that we will need to apply yet 
another transformation, namely T4301 from the list, to the resulting series. Finally, a doubling 
formula for rising factorials, see e.g. [22 1.25], has to be applied. We thus have: 

^ f-k,k + l,-n \ ^ f -n,-k/2,k/2 + l/2,a + n ^ 



l,a J V l,a/2,a/2 + 1/2 

(a/2 + fc/2,a/2 + A;/2 + l/2)„ / -n, -/c/2, 1/2 - /c/2, 1 - n - a 



(a/2,a/2 + l/2)„ ^ V 1> 1/2 - ^ - a/2 - A;/2, 1 - n - a/2 - A;/2 ^ 

{a + k)2n ( -n,-k/2,l/2-k/2,l-n-a 



a)2n Vl>l/2-«-a/2-A;/2,l-n-a/2-A;/2^ 



This shows that formulas (j35|l and (j36l) are indeed equivalent. 

The next step is to compute expressions for {2n\ exp(iX)|2n), by using the series expansion 
for the operator exp(iX). We first use ([35]) : 



(^2nW^\2n) = Y,^-^{MX''\2n) 



i2k)\ 

k>0 ^ ' 



V" M (2i)!(a + n- j)j 2 , 2-,2j f a + 2n + ^2 



j=0 



2i + l/2' 



(37) 



This simply follows from swapping the summation order and writing the inner sum again as a 
hypergeometric series. The same procedure is applicable to formula (I36p . but here one gets: 



Basically, all that is left to do is to integrate these expressions against the complex exponen- 
tial exp(— i(Ap + ^q)). For the expression (j37|) . we proceed as in the case of the ground state. 
We apply Kummer's transformation ([7]) to each of the confluent hypergeometric series. One 
will then get a double summation over integrals of the type (j3ip . (Here, we thus again assume 
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that the order of summation and integration may be switched, which is certainly the case when 
the parameter a is of the form (j30p .) Apphcation of (j3ip for each of the integrals in the series, 
will yield a triple sum expression for W2n(p, q), with one summation "hidden" in the Laguerre 
polynomial. Explicitly, we have the following result: 

As a sanity check, we would like to see that for a = 1/2 formula ([^^ reduces to W2n{p,q)- 
This means that the summation part of (|39p should equal L2n(2p^+2g^). This is not immediately 
clear, so we substitute a = 1/2 and write the Laguerre polynomial as an explicit summation and 
the triple sum then becomes: 

(_2n)fc(-fc);(ff^ + g^y f-k/2,-k/2 + l/2 \ 
^ {1/2),IP ^^^V 1 '7' 

where the summation over j is now performed first. The Gauss hypergeometric series, which 
terminates because either —k/2 or —k/2 + 1/2 is a non positive integer, can be summed in a 
unified way by applying transformation [28], T2140] first: 

Apply this summation, and one gets a double summation. The summation over the index k can 
then be summed using the binomial theorem, and the result follows, showing that for a = 1/2 
expression (|39p reduces correctly to the known expression (j24p . 

Can we also do the integration starting from (j38p ? Clearly, we need a Kummer-like trans- 
formation for 2-^2-series that introduces a quadratic decaying exponential. Fortunately, such a 
transformation was published quite recently by Paris [29]. His main transformation expresses 
one 2-^2-series as an exponential function multiplied by an infinite series containing 2-^2-series. 
We, however, only need a simpler version where a finite sum of 2-^2-series is involved: 

As is stands, expression (p8|) is a double sum. After applying (00]), it will be a triple sum, 
and since integrating using (|3ip introduces an extra summation, the resulting expression will 
be a four-fold summation. Since this computation is quite similar to the ones we have already 
done, we just state the result: 

There now exist quite a number of formulas relating sums of Laguerre polynomials to another 
Laguerre polynomial, one of them being the following [301 4.4.1.5]: 



E(-i)^-(^)4i(-) = (-ir4r^(-)- 

fc=0 ^ ^ 
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This identity is also valid for a = 0, provided one defines the Laguerre polynomials as in ([5]). 
Application of this identity then leads immediately to: 



Also in this case, reduction to the known result for a = 1/2 can be verified, and is easier than 
in the previous case, since the identity needed is recorded as |30l 4.4.1.16]. 

We thus have found closed expressions for W2n{p-,Q) namely ([35|) and dHJ. Prom this we 
now immediately find the corresponding expressions W2n+i{p^Q)- We state the following: 

(42) 

This follows immediately from the fact, shown in Appendix B, that 

(2n + l\X^\2n + 1) = (2n|X'^|2n) U->a+i, (43) 

together with the observation that this is the only place in which the parameter a plays a role 
in the computation of the Wigner distribution function. 

6 Discussion and conclusions 

Parastatistics is one of the alternative approaches in particle statistics proposing that there 
are particles obeying statistics different from Bose-Einstein or Fermi-Dirac [23]. It involves an 
additional parameter (often referred to as the order of statistics) , and for one particular value of 
this parameter one recovers the common Bose-Einstein or Fermi-Dirac statistics. In the current 
case, where we just considered one parabose oscillator, this parameter is a, and for a = 1/2 one 
is reduced to the canonical quantum oscillator. 

The simple case treated here can also be viewed as a deviation from canonical quantum 
mechanics. The one-dimensional quantum oscillator is treated as a Wigner Quantum System, 
and allows several representations of osp(l|2) labelled by a positive number a. In some sense, 
a (or a — 1/2) can be viewed as a quantity measuring the deviation from canonical quantum 
mechanics, see (fT6]l . 

We have now obtained, for the first time, explicit expressions of the Wigner d.f. for the 
parabose oscillator stationary states. This allows one to consider "phase space densities" for a 
"parastatistics quantum system" or for a "non-canonical quantum system". 

Let us now examine some plots of the Wigner distribution functions. First of all, note that 
the Wigner d.f. Wn{p,q) is, just as in the canonical case, a function of + only. So the 
3D-plots of Wn{p,q) in (p, g)-space are rotational invariant. This is why it is more convenient 
to give only two-dimensional plots of Wn{p, q) as a function of r, where 

r = \/jP+q^. 

In Figure 1 we have plotted the ground state Wigner d.f. Wo(p, g) of the one-dimensional 
parabose oscillator for values of the parameter o equal to 1/2, 3/2, 5/2 and 7/2 (or m = 
0,1,2,3). One can see that for a = 1/2 (m = 0, the well-known canonical case) one finds the 
usual Gaussian behaviour of the distribution function. As a (or m) increases, the shape changes 
quite significantly. For values of a larger than 1/2, the shape of Wo{p,q) is rather similar to 
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the shape of the Wigner d.f. for the first excited state of the canonical harmonic oscillator. In 
fact, it follows from that Wo{p,q) for a = 3/2 is equal to Wi{p,q) for a = 1/2, so in that 
specific case it is identical to the canonical Wigner d.f. of the first excited state. As a increases, 
the maximum value of WQ{p,q) moves further away from the origin. The ground state of the 
parabose oscillator behaves as the first excited state of a canonical quantum oscillator but with 
increasing energy as a increases. 

In Figure 2 we plot the distribution functions Wn{p, q) for a fixed a- value (namely a = 3/2), 
and for a number of stationary states (n = 0, 1, 2, 3). Note the separation of even and odd states 
near the origin. 

The properties of the Wigner distribution functions are in agreement with the "wave- 
mechanical representation" for the parabose oscillator, as discussed in [241 Chapter 23]. The 
wave functions in the position representations have also been given here in (jl7l) . Note that these 
wave functions are also the solutions of a wave equation of the form 

The even solutions are the functions (?) with a = g + 1/2 and E = g + \/2 + 2n, and the odd 
solutions are the functions ^2n+i(9) '^ith a = g — 1/2 and E = g + 1/2 + 2n. Equation is 
often referred to as that of the non-relativistic singular harmonic oscillator, or of the Calogero- 
Sutherland model |3H l32 l [33]. The singular behaviour of the last term of the Hamiltonian in (j44p 
is once more in agreement with the negative values of the Wigner distribution function near the 
origin. 

To summarize, we have in this paper constructed the Wigner distribution function for the 
one-dimensional parabose oscillator. This is a quantum system for which the Hamiltonian is 
given by the standard expression ([8]), and for which the position and momentum operators p and 
q do not satisfy the canonical commutation relations, but instead the more general compatibility 
conditions ([9]). This non-canonical situation makes the definition (and computation) of the 
Wigner function more difficult. The approach of Section 3 allowed us to define the right Wigner 
function, satisfying all necessary properties. The computation of the Wigner function Wn{p,q) 
for stationary states is still quite involved, and uses some techniques from hyper geometric series. 
The main results of the computation are given in (j39p . (j4ip and (j42p . On the basis of some 
plots, the behaviour of the Wigner distribution function is further discussed in the context of 
parastatistics or in the context of the wave-mechanical representation of the parabose oscillator. 
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Appendix A 

We have to show that 

i // e-(^'+^')/'e-^(^^'+'^^)(A2 + dXdfi = ^e-P"-'^"k\ 4^ Lk{p^ + q^), (45) 
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Figure 2: The Wigner function of the stationary states of the one-dimensional parabose oscillator 
Wnip, q), for fixed value of a = 3/2 (m = 1) and n = 0, 1, 2, 3. 
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for k G Z+ and p and q real. Let us call the left hand side of the expression ([^5]) /, and apply 
the binomial theorem to the factor (A^ + fj,"^)^: 

Each of the two simple integrals above has a known expression in terms of Hermite polynomi- 
als [34, 2.3.15.10]: 

yielding 

Finally, there is a known summation formula for Hermite polynomials [34^ 4.5.2.4]: 

('')H2,{x)H2k-2j{y) = (-l)'^A;!4^L,(x2 
Application of this summation formula immediately yields the desired result ()45p . 

Appendix B 

In this Appendix, it is our aim to prove formula (136p . We will do so by given formulas for the 
matrix elements {2n + 2l\X'^^\2n) , with I G Z, which will be seen to reduce to formula (j36p when 
I = 0. The proof relies on the fact that the formulas given below satisfy the same recurrence 
and the same boundary conditions as the matrix elements themselves. 
With X as before, one sees that 

= (a+)2(6+)2 + a+Q-{6+, !>-} + {a-f{b-f, 

and hence one gets, using ([HI) and (fT5]l : 

(2n + 2/|X2*=|2n) = (2n + 2l\X'^''-'^ X'^\2n) 

= {a+ f{2n + 2l\X^''-^{b+f\2n) + a+a-(2n + 2l\X'^''^^b+ , b-}\2n) 

+ (Q-)2(2n + 2l\X^^-^{b-f\2n) 
= 2{a+f^/{n + a)(n + 1) (2n + 2l\X^''-^\2n + 2) 
+ 2a+a~(2n + a) (2n + 2/|X2'="2|2n) 

+ 2(a")V"(" + a- 1) (2n + 2/|x2''"^|2n - 2). 

So, if we now denote (2n + 2^|X^'''|2n) as Fk,i{n), the recurrence that Fk^iin) satisfies is: 

Fk,i{n) = 2(a+) V(n + + 1) + 1) + 2a+a-(2n + a) Ffc_i,i(n) 

(4bj 

+ 2(a-) V™(" + a - 1) - 1). 



16 



From the fact that we are deahng with orthonormal states |n), it follows immediately that 



Fofl{n) = 1, (47) 

and also from the actions ()14p and (jlSp it is clear that 

Ffc,;(n) = 0, if |/| > k. (48) 

These two boundary conditions determine the unique solution of the recurrence (j46p . We will 
now give explicit formulas that satisfy the boundary conditions (j47p and (j48p , and we will verify 
that these formulas also satisfy (j46p . The verification of that last fact will boil down to checking 
that two polynomials are equal. As a polynomial basis, we will however not choose the classical 
basis x^, but rather we will work with the polynomials {x)^. 
We now state the formulas, 



(2n + 2/|12fc|2n) = (-l)'2^(a+)'=+'(a-)'^- + + a)i 

(a + / + j)k-i-ji-n)j{-k)i+j{k + I + I) J (49) 



m'm{n,k—l) 

E 

.3=0 



and 

(2n - 2l\X^''\2n) = (-l)'2'=(a+)^-'(a-)^+V(?^ " I + -l + a)i 

min(n_^,fc-0 ^ ^ ^ j)^_i_^{-n + l),{-k)i+j{k + / + 1),- (50) 

where, in both cases, / > 0. It is immediately clear that formulas (j49p and (j50p satisfy the 
boundary conditions (f47l) and (l48]) . Also, the reduction to ([36]) for Z = is readily seen. 

Since the explicit formulas for the matrix elements are different depending on I being positive 
or not, checking that (I49p and (I50p indeed satisfy (I46p will involve three cases: I > 0, I = 
and Z < 0. Since all three are very similar, we will only deal with the case I > here. We will 
abbreviate the summation in (f49l) as Sk,i{n), then substituting the expression (j49]l in (f46l) allows 
simplification of the powers of 2 and as well as the square roots, and leads to the following 
that needs to be checked: 



Sk,i{n) = -Sk-i,i-i{n + 1) + (2n + a)Sk-i,i{n) - n{n + a - l)S'fc_i,;+i(n - 1). 

This is a polynomial identity in the variable —n. Let us denote —n as x. The coefficients of 
the polynomials {x)j are immediately clear for the Ihs. For the right hand side, one can use the 
following easy identities: 

(x = {x)j - j{x)j-i, and x{x)j = {x)j+i - j{x)j. 

This allows to extract the coefficient of {x)j in the right hand side. This coefficient is a rational 
function in the variables a. A;, j and I. After clearing the denominator and simplifying the 
remaining rising factorials, one ends up with a polynomial of a low degree that needs to be 
identically zero, and that is indeed the case. 

So, although the checking is rather tedious, it is at the same time also simple and indeed 
shows that our formulas for the matrix elements are correct. 
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We end this appendix by noting that one indeed has: 

(2n + 1 + 2l\X^\2n + 1) = (2n + 2/|X^|2n)|„^„+i. 

Indeed, it is easily checked that the recurrence satisfied by Gk^i(n) = {2n + 1 + 2l\X^\2n + 1) is 
identical to (06]) but with a replaced by a + 1. Also, the boundary conditions for G are identical 
to (jl7|) and (jlSj) and as these do not involve a the statement is now proved. 



Appendix C 



In this Appendix, we are briefly going to discuss some aspects of the convergence of the main 
expressions (|32]) and ([39]) . The main result is that expressions (f32]) and ([39]) converge absolutely 
when a > 1 for all p and q. In the origin, i.e. (p, q) = (0, 0), the condition a > 1 is also necessary 
for convergence. For the odd states, we thus have absolute convergence for the expressions 
derived from (|32p and (|39p everywhere for a > 0. The only remaining question is thus what 
happens for the even states when [p, q) ^ (0, 0) in the case < a < 1 (and a ^ 1/2). 

First, let us recall the well known fact that if one has two series of positive terms 6^ and 
Cfc such that hk < c^, then convergence of Cfe implies that of bk- Secondly, there are 
some well known tests for convergence of positive series, one of them is Raabe's test which can 
be used if the ratio test fails: let bk be a series of positive terms and let: 



lim k( 

k^oo bk+l 



1) = L, 



one then has that 5^ converges if L > 1 and diverges whenever L < 1. When L = 1, the test 
is inconclusive. 

Let us abbreviate p'^ + q'^ as t. The convergence of ([32p then reduces to the convergence of 



E 

fe>0 



(l/2-a)fc 

:i/2)fc 



(51) 



fc>0 



Next, we are going to use an inequality for the Laguerre polynomials Lk{t) which gives an upper 
bound for \Lk{t)\ not depending on k [35l 10.18(3)]: 



\Lk{t)\ < e*/^ for t > 0. 



(52) 



This gives us: 



(l/2-a)fc 



(1/2)a 



Lkit) 



< 



(l/2-a)fc 



(1/2)a 



The ratio test fails on the series c^, but Raabe's test gives: 
Ck ,^ ,Y |1 + 2A;| 



lim k{- 



1 



lim k 



1 



1 



lim k 

k-^oo 



Ck- 



1 + 2k 



l-2a + 2k 



1 



a. 



2a + 2k\ 

So clearly, the series Ylk '^k converges if a > 1, and hence the series ([5T]) is absolutely convergent 
in that case. Note that the fact that Ylk diverges for a < 1 does not imply anything about the 
convergence of ([32p . The issue about the convergence of (|32p thus remains open for < a < 1. 
When t = 0, however, something more can be said: 



VFo(0,0) 



1 ^ (1/2 - a)k 



vr 



A:>0 



(1/2) 



F ('^2-0,1 
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A Gauss hypergeometric series 

converges if and only if Re(c — a — h) > [22j, which in this case amounts to a > 1. 

The convergence behaviour of Wi{p, q) is now easy to determine from ()34p . It is immediately 
clear that Wi{p,q) is absolutely convergent for a > 0. 

For the convergence of (I39p we follow the same strategy. First, it is clear that the convergence 

of 

is a sufficient condition for the convergence of (I39p since the summation over j is finite. Then, 
we again use ([52]) . and in this case we get 



\bk\ < 



k\ (1/2 -a- 2n)k 



(l/2),(A;-2j)! 



e*/2 = c 



Raabe's test yields in this case that the series Ylk is convergent if a > 1 + 2j — 2n. Since 
j < n (see the summation bounds in ()39p ). a sufficient (and necessary) condition for this to be 
true is a > 1. For t = 0, use of the convergence of the Gauss hypergeometric series for z = 1 
again yields convergence if and only if a > 1. Finally, for the functions W2n+iip, q), we will thus 
have (absolute) convergence everywhere for a > 0. 
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